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Abstract 

We investigate the cosmology of a class of model with noncanonical scalar field and 
matter in an anisotropic time dependent background. Writing the Einstein Equations in 
terms of dimensionless dynamical variables appropriately defined for bouncing solutions, we 
find all the fixed points. From the bouncing conditions and stability of fixed points, solutions 
describing non singular bounce are obtained. 


1 Introduction 

There are two scenarios exist in literature, namely, inflation and bouncing model which address 
the shortcomings of the standard model of cosmology. Though inflation solves most of the 
problems (horizon, flatness and entropy) of the standard model of cosmology, the issue with 
the initial singularity is not resolved under its domain. It is the alternate scenario, nonsingular 
bouncing model, that eradicates the singularity by constructing a universe which begins with a 
contracting phase and then bounces back to an expanding phase through a non zero minimum 
in the scale factor. Nonsingular bouncing models can be categorised into two types, matter 
bounce model [1] and and Ekpyrotik models [2j. For a review on these scenarios refer to [3J 
and [3]. 

A severe problem with bouncing cosmologies is that of instability which develops due to 
the growth of anisotropic stress during contracting phase, as it grows as sixth inverse power 
of scale factor. This problem is known as BKL instability [5]. Originally Ekpyrotic bouncing 
models were developed to cure this instability [6j. This requires a matter field with equation 
of state parameter w greater than 1. In this case energy density of the matter field dominates 
over anisotropic stress so that this instability is eliminated. 

It may be impossible to build a consistent Ekpyrotic bouncing model with one canonical 
scalar field which is free from BKL instability [2] • The Ekpyrotic scenario given in [2] is not able 
to produce scale invariant power spectrum with one scalar field. However, the scale invariant 
spectrum can be achieved by non trivial matching of curvature perturbation across the bounce 
in an Ekpyrotic scenario [7]. Later, a New Ekpyrotic scenario has been developed, wherein 
a second scalar field is added which is an isocurvature mode to start with and becomes an 
adiabatic one during the evolution [8]. This New Ekpyrotic model produces a scale invariant 
power spectrum. BKL instability develops also in New Ekpyrotic models [TO]. It is shown that 
the ’’New Ekpyrotic” scenario can produce nonsingular bounce if the second scalar field is a 
ghost condensate [9j. Unstable growth of curvature perturbation and anisotropy is still a matter 
of concern for New Ekpyrotic models [10] . Whereas matter bounce models can be free from 
these instabilities as shown in ref. Em 

‘email: sukanta@iiserb.ac.in 

^email: manabendra@iiserb.ac.in 


1 



In this paper we consider a noncanonical scalar field with a general function of kinetic term 
F(X), where X = —l/2d ll (j)d fJ '(j). These theories are originally motivated to provide a large ten¬ 
sor to scalar perturbation in inflationary settings H2HH] ■ Dark energy with a general kinetic 
term F(X) is modeled first in ref. 1151 . For other variants of models of dark energy in this 
context refer to m- Other works related to unifying dark matter, dark energy and/or inflation 
for noncanonical scalar field models are studied in HZH203. In order to study the phase space 
in our model, we write the first order equations of motion in terms of dimensionless dynam¬ 
ical variables [21j . The motivation to use noncanonical scalar field as matter is to construct 
nonsingular bouncing models. The phase space analysis of a cosmological model with scalar 
field Lagrangian F(X) — V((j >) and matter for an FRW background is given in ref. [221 . The 
condition for nonsingular bounce is also discussed in ref. [[22]. In order to explore the growth 
of anisotropy(shear parameter) near bounce and late time isotropization, we do a phase space 
analysis for the model discussed in [ 22] with a Bianchi I metric. This can be easily extended 
to other nonsingular bouncing models. This analysis may help us to avoid BKL instabilities 
in bouncing models. It is shown that growth of initial anisotropy and inhomogeneity can be 
chosen at a level consistent with observations which will be useful for checking the viability of 
bouncing cosmological models [23]. 

Here, we study the cosmology of an anisotropic universe with a matter Lagrangian of the 
form F(X) — V(<j>) and an additional matter. In section [2] we write the Einstein equations in 
terms of dynamical variables suitable for studying bouncing scenarios. Fixed points and their 
stability are analysed in section [3] Conditions for existence of nonsingular bouncing solution, 
in terms of dynamical variables are derived in section [4j We summarise our results in section [5j 

2 Einstein Equations in Bianchi I background 

The action for our model is given by 

S = J d^x^gi-.R + F(X) - V(<j>) + L m \ ( 1 ) 

where L m , R, g and <f> represent the Lagrangian of the matter field, the Ricci scalar, determinant 
of the metric and the scalar field respectively. 

To see the behaviour of anisotropy of the spacetime in a nonsingular bouncing scenario we 
choose to work with a homogeneous but anisotropic metric(Bianchi I) with a planar symmetry. 
The line element of Bianchi I metric with planar symmetry is 

ds 2 = — dt 2 + a 2 (t)dx 2 + b 2 (t)(dy 2 + dz 2 ). (2) 

We define average hubble parameter H and Shear h as 


H = \{Ha + c 2H h ) 

, H b - H a 

h = --=— 

\/3 


(3) 


In terms of averaged Hubble parameter and shear h, the Einsteins equation take the following 
form 


d H 

dt 
d h 

dt 

H 2 


—H 2 — ^h 2 - p(/> + 3p), 
3 fa 

-3 hH , 

p h 2 

-A - 

3 ^ 3 ’ 


(4) 
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where p = p ( p + p m and p = P(f) + p m . 

Here the energy density p 0 and pressure p<f> of the scalar field is found to be 


Pp = 2XF x — F + V , 

P* = F(X) — V (</>), (5) 

and p rn and p rn are the energy density and pressure due to L m . 

Substituting Eq. © in first and third line of Eq. Q , we get 


^ = -H 2 + \h 2 - \{2XF x - F + V + p m + 3(F — V) + 3p m ) 
dr 3 o 

and 

r2 = 2 XFx - F V Prn 

3 ^ 3 ^ 3 ^ 3 


( 6 ) 

(7) 


Here we further define few more variables which are useful for defining dimensionless dy¬ 
namical variables. They are 


Pk = 2 XFx — F, 

F 

Wk ~ 2 XF x - F ’ 

1 d logV 

V?>\pk\ dt 


( 8 ) 


where p k is the kinetic part of the energy density p^, w k is the ratio of kinetic part of the 
pressure to the p k and a is the auxiliary variable which depends on the variation of potential 
with time. 

Neglecting the interaction between scalar field and matter, continuity equation for p^ in 
terms of dimensionless time variable N (dN = Hdt ), is 


A^(2XF x -F + V) + 6 XF x =0. (9) 

Now we define a set of dimensionless dynamical variables which is suitable for nonsingular 
bounce models. Relevance of these variables are that they remain finite during the entire 
evolution across bounce. The dynamical variables are 


x = 


\/3 H_ 

vW’ 


y = 


VW\ 


signiV ), z = 


Using Eqs.flTUJ), ([!]) and 
y and 5 are written as, 


vW 

and parameters defined in Eq. 


7 ’ q _ Pm 

vW’ m “l pk\' 


( 10 ) 

, the evolution equations of x, 


dx 

dfV 

d y 

dJV 

d5 

dN 

dH m 

dN 


= ~2 “ Wm)sign(pk) + (1 + w m )(x 2 - y\y\) + (1 - w m )z 2 

3 

+ ^x[(w k + l)x - ay\y\sign( P k)], 

3 

= + (w k + l)x - cry\y\sign{ P k)\, 

= -3 zx + 3zx(l + w k ) - 3zy\y\sign(p k ), 

= -3(1 + w m )xQ m - [3ay\y\sign(p k ) - 3x(l + w k )\ , 


( 11 ) 
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where dN = y and the constraint equation relating dynamical variables is 

x 2 ~ y\y\ ~ z 2 ~ &m = 1 x sign(pk). (12) 

The equation for parameter a becomes [22] 

iE = -3cr 2 (r- 1) + 3ct (2S (Wk + 1 ) + Wk ~ l) 
dN 2 ( 2(7 + 1 ) (w k + 1 ) 


{w k + 1 ) x - ay* 


(13) 


where S = X p xx and T = V y^ . 

For our model we have taken power law form for F(X) = FqX v , where Fq is a constant. 
For this form of F(X), w k = aR d S = rj — 1. 

Potential V(4>) is taken as V (</>) = Voe _c ^, where Vo and c are constants with positive values. 
For this choice of V (</>), T becomes unity. 

In the next section, we do a fixed point analysis of dynamical equations for x, y, z and a. 
The evolution of Q m is determined from the constraint Ea. ([12l) . 


3 Fixed Point Analysis 

In this section, we do a fixed point analysis of our system of dynamical equation in order to 
extract the qualitative information about the nature of solution. Fixed points are calculated by 
taking the first derivative of the dynamical variables to be zero. The stability of a fixed point 
is determined from the behaviour of a small perturbation around that fixed point. 

We get the set of fixed points x c , y c , z c and a c by solving the following set of equations 
simultaneously (where the subscript c denotes fixed points). Now, if we define the slopes of the 
dynamical variables x, y, z and a as f(x,y,z,a), g(x,y,z,a ), h(x,y, z,a) and i(x,y,z,a). The 
set of equations we need to solve to obtain the fixed point is 


where, 

f(x,y,z,cr) 

g(x,y,z,a) 
h(x,y,z,a) 
i(x,y,z,cr) 


f(x,y,z,a) 
g(x,y,z,a) 
h(x, y, z, a) 
i(x,y,z,a) 


4i = 0 , 

dIV 

«- 0 , 
d N 
dz 

= 0 , 
dN 
da 

= 0 , 

dN 


(14) 


-^[(wfc - w m )(signp k ) + (1 + w m )(x 2 - y\y 

3 

+ ( Wk + 1 ) 5 “ ^y\y\{signp k )\ , 

-3 zx + 3zx{l + w k ) - 3zy\y\sign(p k ), 

3 [2z>(w k + 1) + (w k — 1)] u , i\~ ~2i 

2 2(to+l)K + l) 1(wt + 1)x - ay t 




The corresponding fixed point for Cl m can be found using the constraint Eci. (H2l) . 

The stability of the fixed points can be examined from the evolution of small pertubations 
around fixed points. Now, if (x c , y c , z c , a c ) is a fixed point and Sx = x — x c ,5y = y — y c , 


ay\y\sign{p k ) 


(1 
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5z = z — z c and 5a = a — a c be the respective perturbation around it, then the evolution of the 
perturbation is determined by 


5x = x = f(x c + 5x,y c + 5y,z c + 5z,a + 6a ), 
5y = y = g(x c +6x,y c +Sy,z c +5z,a + 5a), 
6z = z = h(x c + 6x, y c + 5y, z c + 6z, a + 6a ), 
6a = a = h(x c + 6x, y c + 6y, z c + 6z, a + 5a ), 


(16) 


The evolution equations,upto first order, for these pertubations are 


( 5x ^ 


( 5x \ 

5y 

= A 

5y 

5z 


5z 

\6& ) 


\5a ) 


where the matrix is 


A = 


/ ^ 

dx 

dg_ 
dx 
dh 
dx 
di 
dx 


9f 

dij 

9g 

dij 

dh 


1 

dy 


d£ 

dz 

dg_ 

dz 

dh 

dz 

di 

dz 


9£ 

da 

dg_ 

da 

dh 

da 

di 

da 


\ 

/ 


(17) 


(18) 


is the Jacobian matrix and is evaluated at the fixed point (x c , y c , z c , a c ) and hence each entry 
of A is a number. The solution of the system of equations can be found by diagonalizing the 
matrix A. A non trivial solution exists only when the determinant |A — AI| is zero. Thus, 
solving this equation in A we would get all the eigen values of the system corresponding to each 
fixed points. 

We have two cases: one with positive kinetic term, sign(pk) = +ve and other one with 
negative kinetic term, sign(pk) = —ve.. 


3.1 Case I, sign(pk) = +ve 


In this case, we study the fixed points for all possible values of parameters. The fixed point 
(0, 0,0, 0) is obtained for Wk = w m signifying all the dynamical variables x, y, z and a, going 
to zero at late times. It is a nonhyperbolic fixed point as the eigen value of A for this is 
(0,0, 0,0). It’s stability cannot be decided from our first order analysis of perturbations. From 
now onwards eigenvalues would mean eigenvalues of matrix A for the rest of the paper. 

The second fixed point (1,0,0, 0) denotes a late time kinetic dominated universe with other 
dynamical variables y, z and a becoming zero. In this case, eigenvalue is ( 3 ( w ^ +1 ) , —3 + 3(1 + 
w k), f (—1 + Wk+ ( ' 1 ~ 1 J ' fc ^ 1 +Wfc ' > ))3(l + u;fc) — 3(l + u' m )). This is a stable fixed point for the region 
of parameter space shown in the Fig. jT]. 

The next stable fixed point is (—1,0,0,0) with eigenvalue (|(—1 — Wk), 3 — 3(1 + Wk), 


3 (_l_ Wfc )(_l +1 ^ + ( 1 ~ m k)( 1 + w k) ) 

- 2 ( 1 +Wk ) - — -> f (—1 — Wk) — §(1 + Wk) + 3(1 + w m )) shows again a late time 

kinetic dominated phase but with a negative value of averaged Hubble parameter H signifying 
a contracting universe. This fixed point is found to be stable for the region of parameter space 
shown in Fig. [2]. The point (—1,0,0,0) may not be important for bouncing point of view, as 
we need the universe to transit to an expanding phase to be discussed in section 01 

The next three fixed points being (0, 0, —^4=^, 0) for Wk > w m and w m > 1, (0, 0, , 0) 

V L-\-Wrn V 

for w k > w m and w m > 1 , ( 0 , ^^^^, 0 , 0 ) with eigen values ( 0 , 0 , -3 - w k w m , 3 - w k w m ), 
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Fixed Points (x c , y c , z c , a c ) 

Stability Conditions 

(0,0,0,0) for w k = w m 
(1,0,0,0) 

(-1,0,0,0) 

(o, o, J_\ +W 2 ’ °) Wlth w k > w m and w m > 1 
(0, 0, j_\ +w 2, 0) with w k > w m and w m > 1 

0. .0.0) 

Can’t decide 

Stable for w k < —1 and w rn > 0, see Fig.pQ 
Stable for w k < — 1 and w m > 0, see Fig. [2] 
Can’t decide 

Can’t decide 

Can’t decide 


Table 1: Stablity Analysis of fixed points for sign(pk) = +1 



-10 


10000 


Figure 1: allowed region of parameter space for the fixed point (1,0,0,0) 


(°’ 0 ’-ySSM)’ (°’ -3yjw k + W 2 k -w m - w k w m y/2,3 Wk+u ’ k ^ WkWm ) are 

hyperbolic points. The stability of such fixed points goes beyond the linear stabilty 
All the fixed points and their stability are noted in tablepQ. 


also non¬ 
analysis. 


3.2 Case II, sign(p k ) = —1 

In this section, we state the results of stability anlysis of our dynamical variables for the negative 
sign of kinetic energy density. The fixed points are found to be (0, 0,0, 0), (0, 0, — 0), 

(o, o, , 0) and (0, 0) with eigen values (0, 0,0, 0), (0, 0, -3y/-w 2 k + w k w m , 

3^/-w 2 k + w k w m ), (0, 0, —3a/— u ,2 fc + w k w m , 3a J-w 2 k + w k w m ) and (0, , 0,0) respec- 

tively. All these fixed points are nonhyperbolic and tabulated in table [2]. 


4 Bouncing Scenario 

Now, we obtain the conditions for nonsingular bounce to occur and also show the evolution 
of dynamical variables numerically. A nonsingular bounce is attained whenever the universe 
passes from a contracting phase to an expanding phase through a minimum value of the avearge 
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Figure 2: allowed region of parameter space for the fixed point (—1,0,0,0) 


Fixed Points (x c , y c , z c , cr c ) 

Stability Conditions 

(0, 0,0, 0) for w k = w m 

10 0 ~ w k+Wm Q\ 

fu,u, _ 1+Wm , Uj 

f0 0 - W k+Wm Q\ 

10,0, _ 1+Wm ) U J 

fO ~ W k+Wm Q Q\ 

1 U ’ L+w™ ’ u ’ u f 

Can’t decide 

Can’t decide 

Can’t decide 

Can’t decide 


Table 2: Stablity Analysis of fixed points for sign(p k ) = — 1 
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Figure 3: Evolution of the dynamical variables x (top left), y (top right), z (bottom left) and 
a (bottom right) for the fixed point (x c ,y c ,z c ,a c ) = (1,0,0,0) with the values of parameters 
Wk = —2.0, w rn = 1/3 and sign(pk) = +ve for different initial conditions 







scale factor A(t)(= (aft 2 ) 1 / 3 ) but not zero. Mathematically, it satisfies 

where subscript b denotes value of the variable at the bounce, and 


( 19 ) 




( 20 ) 


for minimum to occur. This implies 


fdH\ _ [A 


A 


V d t J i 


A 


> 0 


( 21 ) 


Now, writing the above conditions in terms of dynamical variables for bouncing, we get 
Xb = 0 and > 0 which translate to the following equations 


/ dx \ { 

\dNj b ~ 


2 L 


and 


(w k - w m )(signp k ) + (1 + w m )(-y\y\) + (1 - w m )z 2 

y\y\ z 2 


> 0 


(1 - w m ) (1 + w m ) 


> 1 x sign(p k )( 


U>k ~ W r 
1 — 1C 2 


Thus at the bounce, we obtain the constraint equation among dynamical variables as 
(x 2 - y\y\ - z 2 - &m) b = ~y\y\ - z 2 - O m = 1 x sign(p k ). 


( 22 ) 


(23) 


(24) 


For different negative initial conditions of x (contracting phase), Fig. [3] ( top left) shows 
its transition to positive values (expanding phase) crossing zero (bounce). The bouncing is 
guaranteed by the positivity of the slope of x as shown in Fig. [5] (left). Fig. [3] ( top left) and 
Fig. [5] (left) do indeed represent a stable bouncing scenario. This is obtained by setting the 
values of equation of state parameters w k = —2, (77 = 1/4), w rn = 1/3 and sign(pk) = +ve and 
sign(y) = +ve. The evolution of other dynamical variables can be noted from Fig. [3], which 
show their asymptotic evolution to the respective fixed points. 

It can be seen that the fixed point (x c ,y c ,z c ,a c ) = (1,0,0,0) does give rise to a stable 
bouncing universe provided it satisfies Eas. fEHl [231) • From this analysis, we conclude that, 
finally, after the bounce our universe is driven by kinetic energy density at late time. The 
nonsingular bounce happens only for negative values of Cl m with our choice of parameters as 
shown in Fig. [ 6 ]. The other fixed point (—1,0, 0,0), though stable, can not give rise to a bouncing 
scenario as it ends up with a negative value of averaged Hubble parameter, H, signifying a late 
time contracting phase. 

One main point of our work is to show the behaviour of shear parameter, z, in a nonsingular 
bouncing set up. The shear parameter increases initially as the universe contracts and then 
decreases to zero value in the expanding phase after crossing the bounce Fig. [3] . Shear param¬ 
eter remains finite at bounce. Hence, no anisotropic collapse as expected in nonsingular bounce 
models. Shear may not dominate near the bounce and BKL instability can be avoided which 
in turn constrain the allowed initial values of dynamical variables. A late time isotropization of 
the universe is attained in our model which is useful for building realistic model. 

Finally, we show the effect of 77 parameter on the behaviour of bouncing solutions in Fig. 
[5] (left). All the plots are generated for the same set of initial conditions and set of parameters 
w m = 1/3, sign(pk) = +ve but with three different values of parameters 77 =1/4,1/6 and 1/8. 
The values of 77 lies between 0 and 1/2 which corresponds to w k < — 1 allowed by the stability of 
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Figure 4: a e vs i on (left) and 5 vs x on (right) for Wk = —2.0, w m = 1/3 and sign(pk) = 
+ve 




Figure 5: x vs N on (left) and z vs x on [right) for 77 = 1 / 4 , 1/6 and 1/8 with w rn = 1/3 and 
sign(pk) = +ve 

fixed points. Thus, higher powers of kinetic term p more than or equal to 1/2 are ruled out by 
stability criteria. It has been observed that the value of p has a direct impact on the occurence 
of bouncing point. Indeed, the position of bouncing point is delayed as we decrease the value 
of p. The right hand side of Fig. [5] shows the variation of shear parameter z with respect to x 
for the three aforesaid values of p. Though, the universe isotropizes at late time in each of the 
three cases, it is seen that the value of shear at the bounce is more for higher values of p. 

5 Conclusion 

A cosmological scenario with a noncanonical scalar field and matter is explored in this work. 
Using evolution equations for a set of dimensionless dynamical variables, we find all fixed points 
for both positive and negative kinetic energy density terms. Stability conditions are obtained 
and allowed region of parameter space is shown. We evolve dynamical variables for different 
initial conditions with a choice of values of parameters allowed by stability conditions. The 
necessary and sufficient conditions for a nonsingular bounce to occur in terms of the dynamical 
variables are derived. Cosmological solutions satisfying both nonsingular bouncing conditions 
and stability criteria are obtained for the choice of parameters. This is achieved for the negative 
energy density of matter with equation of state parameter w m = 1/3. In addition to this, the 
shear parameter is found to be finite at the bounce as expected in a nonsingular bouncing 
scenario. Late time isotropization occurs in the expanding phase after the bounce which is 
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Figure 6: Cl m vs N for Wf. = —2.0, w m = 1/3 and sign(pk) = +ve 


independent of initial conditions. Finally, the effect of the parameter on the behaviour of 
bouncing solution is noted. It is seen that the point of occurence of bounce is delayed as we 
decrease the value of g. 

We restrict our analysis only to positive sign of potential. It is straightforward to extend 
our analysis for a negative potential by changing the parameter sign(y) to be —1. It would 
be interesting to see the growth of an anisotropic perturbation around Bianchi-I metric in a 
nonsingular bouncing model and it’s signature in cosmological observations. 
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